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The discovery of a novel coronavirus (CoV) as the causative agent of severe acute respiratory syndrome 
(SARS) has highlighted the need for a better understanding of CoV replication. The replication of SARS-CoV 
is highly dependent on host cell factors. However, relatively little is known about the cellular proteome changes 
that occur during SARS-CoV replication. Recently, we developed a cell line expressing a SARS-CoV sub- 
genomic replicon and used it to screen inhibitors of SARS-CoV replication. To identify host proteins important 
for SARS-CoV RNA replication, the protein profiles of the SARS-CoV replicon cells and parental BHK21 cells 
were compared using a quantitative proteomic strategy termed “stable-isotope labeling by amino acids in cell 
culture—mass spectrometry” (SILAC-MS). Our results revealed that, among the 1,081 host proteins quantified 
in both forward and reverse SILAC measurements, 74 had significantly altered levels of expression. Of these, 
significantly upregulated BCL2-associated athanogene 3 (BAG3) was selected for further functional studies. 
BAG3 is involved in a wide variety of cellular processes, including cell survival, cellular stress response, 
proliferation, migration, and apoptosis. Our results show that inhibition of BAG3 expression by RNA inter- 
ference led to significant suppression of SARS-CoV replication, suggesting the possibility that upregulation of 
BAG3 may be part of the machinery that SARS-CoV relies on for replication. By correlating the proteomic data 
with these functional studies, the findings of this study provide important information for understanding 


SARS-CoV replication. 


The outcome of a viral infection is regulated in part by the 
complex coordination of viral and host interactions that com- 
pete for the control and optimization of virus replication. Vi- 
rus-host interactions are crucial determinants of virus host 
range, replication, and pathology. Studies of virus-host inter- 
actions have advanced understanding of viral and cellular func- 
tion and can provide targets for antiviral development. One 
area in which the importance of host factors is increasingly 
emerging is the replication of positive-strand RNA [(+) RNA] 
viruses. (+) RNA viruses are the largest genetic class of viruses 
and include significant human pathogens such as severe acute 
respiratory syndrome coronavirus (SARS-CoV), hepatitis C 
virus, and West Nile virus. Defining the host factors that gov- 
ern the replication of (+) RNA viruses will enhance our gen- 
eral understanding of their molecular biology and may have 
important implications for the development of novel antiviral 
control strategies. Whereas recent studies show that host fac- 
tors are critical for (+) RNA virus genome replication and 
mRNA synthesis and are targeted by (+) RNA viruses to 
modulate host gene expression and defenses (1, 33), identify- 
ing such factors remains difficult. 

In 2003, a novel coronavirus, SARS-CoV, emerged from 
zoonotic pools of virus in China to cause a global outbreak of 
SARS (10, 31). The SARS-CoV genome encompasses 29,727 
nucleotides, and the genome organization is similar to that of 
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other coronaviruses. The genome is predicted to contain 14 
functional open reading frames (ORFs) (41, 60). SARS-CoV 
genome translation yields two large replicase polyproteins 
(ppla and pplab) that are autoproteolytically cleaved into 16 
nonstructural proteins (nsp1 to -16) by proteases residing in 
nsp3 and nsp5 (21, 23, 53). These 16 SARS-CoV nsps include 
RNA-binding protein (nsp9), RNA-dependent RNA polymer- 
ase (nsp12), helicase (nsp13), RNA synthesis proteins (nsp8 
and nsp14), and several transmembrane proteins (nsp3, nsp4, 
and nsp6), etc. (69). They are the primary constituents of the 
replication/transcription complex (RTC), which is believed to 
be associated with characteristic double membrane vesicles 
(DMVs) derived from modified host cell membranes (57, 62). 
RNA replication is believed to occur on DMVs and uses host 
proteins as part of their replication strategies (57). Thus, iden- 
tifying such host factors and their contributions has long been 
recognized as an important frontier. 

Recent advances in molecular profiling technologies have 
allowed for advances in our understanding of the mechanisms 
of cellular responses to the SARS-CoV infection. Analysis of 
gene expression profiles during viral infection is one of the 
powerful approaches to probe potential cellular genes involved 
in viral infection and pathogenesis (37, 59), but ultimately 
protein expression and posttranslational modification (PTM) 
determine virus replication. Thus, the molecular analysis of 
viral infection would greatly benefit from a proteomics ap- 
proach that combines the advantages of high-throughput anal- 
ysis and the focus on protein levels and modifications. Pro- 
teomic techniques as a powerful research tool have recently 
become available for large-scale protein analysis, and stable- 
isotope labeling by amino acids in cell culture (SILAC) is one 
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of the most effective methods for the simultaneous detection of 
diverse changes in protein expression (49). With SILAC, the 
entire proteome of a given cell population is metabolically 
labeled by heavy, nonradioactive isotopic variants of amino 
acids, thus making it distinguishable by mass spectrometry 
(MS) analysis. Thereafter, two or more distinctly SILAC- 
labeled cell populations can be mixed and analyzed in one MS 
experiment, which allows accurate quantitation of proteins 
from the different cellular states (49). Because of its properties 
of being simple, inexpensive, and accurate, SILAC is being 
used in the life sciences more and more extensively (4, 9, 30, 
58). To date, a small but increasing number of studies have 
used proteomics approaches to investigate various aspects con- 
cerning the infection and pathogenesis of the SARS-CoV and 
the virus-host interactions (24, 67). In spite of these data, it is 
important to establish a comprehensive catalogue of the cel- 
lular factors interacting with the virus RNA that may regulate 
the SARS-CoV replication. 

In our previous work, we developed the first SARS-CoV- 
derived replicon cell line (19). This SARS-associated replicon 
cell line is based on the use of SARS-CoV replicon cDNAs 
generated by reverse genetic techniques. The viral envelope- 
protein coding S, E, and M genes were replaced by a selectable 
marker and a green fluorescent protein (GFP) reporter gene. 
The nucleocapsid (N) gene was retained for efficient replica- 
tion of replicon RNA. Due to the absence of some structural 
genes, no infectious viral particle could be produced by the 
cells. However, since all trans- and cis-acting components re- 
quired for viral RNA synthesis are retained, these partial viral 
RNAs could replicate autonomously in the cells. Our data 
indicate that this replicon cell line can be applied to high- 
throughput screening for anti-SARS drugs without the need to 
grow infectious SARS-CoV (20). 

To generate additional insights into the molecular events 
controlling SARS-CoV replication, in this study, the SILAC 
method was employed to compare the protein profiles of the 
SARS-CoV-derived replicon cells to those of the parental 
BHK21 cells. Many interesting differentially expressed pro- 
teins were identified that potentially play functional roles dur- 
ing virus replication. Further functional studies demonstrated 
that BAG3 plays an important role in SARS-CoV replication. 
By correlating the proteomic data with these functional stud- 
ies, the current results not only provide insights into the mech- 
anism underlying the virus-host interactions but also have di- 
rect implications for drug development for SARS-CoV. 


MATERIALS AND METHODS 


Cells and viruses. The BHK21 (baby hamster kidney) and Vero E6 (African 
green monkey kidney) cell lines were purchased from the American Type Cul- 
ture Collection (ATCC) and maintained in Dulbecco’s modified Eagle’s medium 
(DMEM; Invitrogen) supplemented with 10% fetal bovine serum (FBS; Invitro- 
gen) at 37°C in 5% CO . The SARS-CoV replicon cell line SCR-1 stably ex- 
pressing a SARS-CoV replicon encoding GFP has been described previously (19, 
20), and the cells were cultured in DMEM (Gibco) containing 10% fetal calf 
serum (Gibco) and 10 yg/ml blasticidin (Invitrogen). A seed stock of SARS-CoV 
(strain SIN 2774) passaged in Vero E6 cells was used for infection. The Beau- 
dette strain of the coronavirus infectious bronchitis virus (IBV) (ATCC VR-22) 
was obtained from the ATCC and was adapted to Vero E6 cells as described 
previously (40). 

Protein analysis by SILAC labeling. To differentially label the SARS-CoV 
replicon cell line SCR-1 and parental BHK21 cells, the SILAC protein quanti- 
tation kit (Pierce Biotechnology, Rockford, IL) was used according to the man- 
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FIG. 1. Schematic showing application of forward and reverse 
SILAC combined with LC-MS/MS for the comparative analysis of 
protein expression in SARS-CoV replicon cells (SCR-1) and their 
parental BHK21 cells. 
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ufacturer’s instructions. In brief, cells were grown in SILAC DMEM (Pierce 
Biotechnology, Rockford, IL) containing 10% (vol/vol) dialyzed FBS, and either 
0.1 mg/ml heavy L-[‘*C,]- or Light L-['?C,]lysine (Pierce Biotechnology, Rock- 
ford, IL). Cells were propagated in SILAC medium for >6 generations to ensure 
nearly 100% incorporation of labeled amino acids. In the forward SILAC ex- 
periment, the SCR-1 cells were cultured in light medium, whereas the BHK21 
cells were cultured in heavy medium. Reverse SILAC experiments were also 
performed in which the SCR-1 and BHK21 cells were cultured in the heavy and 
light medium, respectively (Fig. 1). For protein extraction, cells were lysed in 8 
M urea and 20 mM HEPES supplemented with complete protease inhibitor 
cocktail tablets (Roche, Nutley, NJ). All chemicals were purchased from Sigma- 
Aldrich (St. Louis, MO) unless otherwise stated. Cellular debris was removed by 
centrifugation for 30 min at 13,200 x g and at 4°C. Protein concentrations were 
measured in duplicate using the RC DC protein assay (Bio-Rad, Hercules, CA). 
Protein separation and in-gel digestion. The light and heavy cell lysates were 
combined at 1:1 ratio (wt/wt) (100 wg in total), boiled in SDS-PAGE sample 
buffer, separated by 12% SDS-PAGE, and stained with Coomassie blue. The 
entire gel lane was cut into 30 sections for in-gel tryptic digestion. The excised 
sections were chopped into small particles and washed in water and then com- 
pletely destained using 100 mM ammonium bicarbonate in 50% acetonitrile 
(ACN). A reduction step was performed by addition of 100 pl of 10 mM 
dithiothreitol (DTT) at 37°C for 1 h. The proteins were alkylated by adding 100 
pl of 50 mM iodoacetamide and allowed to react in the dark at 20°C for 30 min. 
Gel sections were first washed in water, and then acetonitrile, and finally dried 
with a SpeedVac centrifuge (Thermo Fisher Scientific, Waltham, MA). Diges- 
tion was carried out using 20 pg/ml sequencing-grade modified trypsin (Pro- 
mega) in 50 mM ammonium bicarbonate. Sufficient trypsin solution was added to 
swell the gel pieces, which were kept at 4°C for 45 min and then incubated at 
37°C overnight, after which peptides were extracted from gels with 5% acetic acid 
in H,O and in CH3CN/H,0O (1:1 [vol/vol]). The resulting peptide mixtures were 
dried using vacuum centrifugation and stored at —80°C for further analysis. 
Protein identification and quantification. Online liquid chromatography-tan- 
dem MS (LC-MS/MS) analysis was performed on an Agilent 6510 Q-TOF system 
coupled with an Agilent high-performance liquid chromatography (HPLC)-Chip 
Cube MS interface (Agilent Technologies, Santa Clara, CA). The sample injec- 
tion, enrichment, desalting, and HPLC separation were carried out automatically 
on the Agilent HPLC chip with an integrated trapping column (160 nl) and a 
separation column (Zorbax 300SB-C,; 75 wm by 150 mm, 5-m particle size). 
The peptide mixture was first loaded onto the trapping column with a solvent 
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mixture of 0.1% formic acid in CH3;CN/H,O (2:98 [vol/vol]) at a flow rate of 4 
l/min. The peptides were then separated with a 90-min linear gradient of 5 to 
60% acetonitrile in 0.1% formic acid and at a flow rate of 300 nl/min. The chip 
spray voltage (VCap) was set as 1,950 V and varied depending on chip condi- 
tions. The temperature and flow rate of the drying gas were set at 325°C and 4 
liters/min, respectively. Nitrogen was used as the collision gas, and the collision 
energy followed an equation with a slope of 3 V/100 Da and an offset of 2.5 V. 
MS/MS experiments were carried out in the data-dependent scan mode with a 
maximum of five MS/MS scans following each MS scan. The m/z ranges for MS 
and MS/MS were 300 to 2,000 and 60 to 2,000, and the acquisition rates were 6 
and 3 spectra/s, respectively. Agilent MassHunter workstation software (version 
B.01.03) was used to extract the MS and MS/MS data. The data were converted 
to m/z data files with MassHunter Qualitative Analysis. Mascot Server 2.2 (Ma- 
trix Science, London, United Kingdom) was used for protein identification by 
searching the m/z data files against the IPI mouse protein database (version 3.21; 
51,432 sequences) (28). The maximum number of miscleavages for trypsin was 
set as 1 per peptide. Cys (+57.0215 Da; carbamidomethylation) was set as fixed 
modification, whereas Met (+15.9949 Da; oxidation) and Lys (+6.0201 Da; 
SILAC heavy amino acid) were considered as variable modifications. The mass 
tolerances for MS and MS/MS were 50 ppm and 0.6 Da, respectively. Peptides 
identified with individual scores at or above the Mascot-assigned homology score 
(P < 0.01 and individual peptide score of >30) were considered as specific 
peptide sequences. The false discovery rates (FDR) determined by decoy data- 
base search were <0.95%. All identified peptides were subjected to relative 
quantification analysis using the program Census (51). Only proteins with a 
minimum of 2 quantifiable peptides were included in our final data set. The 
protein ratios were calculated from the average of all quantified peptides. The 
quantification was based on four independent SILAC and LC-MS/MS experi- 
ments, which included two forward and two reverse SILAC labelings, and the 
proteins reported here could be quantified in both forward and reverse SILAC 
experiments. 

Protein categorization. Differentially expressed proteins (DEPs) were classi- 
fied based on the PANTHER (protein analysis through evolutionary relation- 
ships) system (http:/Avww.pantherdb.org), which is a unique resource that clas- 
sifies genes and proteins by their functions (43). 

Western blot analysis. Protein extracts (30 4g) prepared with radioimmuno- 
precipitation assay (RIPA) lysis buffer (50 mM Tris-HCl, 150 mM NaCl, 0.1% 
SDS, 1% NP-40, 0.5% sodium deoxycholate, 1 mM phenylmethyl sulfonyl fluo- 
ride [PMSF], 100 mM leupeptin, and 2 mg/ml aprotinin, pH 8.0) were resolved 
with a 10% SDS-PAGE gel and transferred onto polyvinylidene difluoride 
(PVDF) membranes (Millipore Corporation, Billerica, MA) by electroblotting 
and then blocked using Tris-buffered saline-Tween 20 (TBST) buffer containing 
5% nonfat milk. The membranes were probed with rabbit anti-heat shock protein 
90 (anti-HSP90) polyclonal antibody, goat anti-ENO1 polyclonal antibody, goat 
anti- YWHAZ polyclonal antibody, goat anti-RPS19 polyclonal antibody, rabbit anti- 
HSPAIA polyclonal antibody (Santa Cruz Biotechnology, Santa Cruz, CA), rabbit 
anti-BAG3 polyclonal antibody, rabbit anti-hnRNP A1 polyclonal antibody (Abcam, 
Inc., Cambridge, MA), mouse anti-lactate dehydrogenase B (anti-LDHB) monoclo- 
nal antibody, rabbit anti-cell division control protein 42 (anti-CDC42) polyclonal 
antibody (Abnova, Walnut, CA), mouse anti-tubulin monoclonal antibody (Lab- 
Vision, Fremont, CA), and anti-nsp5 (3CLpro) antibody (20), respectively. After 
washing, the membranes were incubated with horseradish peroxidase (HRP)- 
conjugated secondary antibodies at room temperature for 1 h. The membranes 
were washed thoroughly with PBST for three times. The secondary antibody was 
detected by using the ChemiGlow chemiluminescence reagents (Alpha Innotech, 
San Leandra, CA). Finally the immunoblots were scanned, and densitometric 
analysis was performed using the public domain NIH Image program ImageJ 
(developed at the U.S. National Institutes of Health and available on the Inter- 
net at http://rsb.info.nih.gov/nih-image/). 

siRNA for BAG3 and transfection of siRNA. SCR-1 cells were transfected with 
50 nM small interfering RNA (siRNA) specific to mouse Bag3 or with nontar- 
geting siRNA (Dharmacon, Lafayette, CO). siRNAs were introduced to cells 
using Lipofectamine Plus in Opti- MEM medium (Invitrogen, Carlsbad, CA). 
siRNA has been selected among four different siRNA sequences that have been 
evaluated for high specificity and lack of off-target effect at the concentration 
used. Forty-eight hours posttransfection, fluorescence microscopy was used to 
observe the changes of green fluorescence in SCR-1 cells, protein extracts were 
analyzed by Western blotting, and total RNA was purified and used for real-time 
PCR analysis. 

Establishment of stable cell line expressing BAG3 shRNA. SMARTvector 2.0 
lentiviral short hairpin RNA (shRNA) particles specific to human Bag3 or 
nontargeting control particles (Dharmacon, Lafayette, CO) were transfected into 
Vero E6 cells according to the manufacturer’s instructions. Transfected cells 
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were selected by using 1 g/ml of puromycin. One clone was selected from BAG3 
shRNA-transfected Vero cells (designated as Vero-KD, for “Vero knockdown”), 
and one clone was selected from control shRNA-transfected Vero cells (desig- 
nated as Vero-NC). BAG3 protein knockdown was assessed by Western blotting 
analysis. 

Plaque reduction assay. The plaque reduction assay followed the procedures 
previously described (20, 39). All procedures involving manipulation of live 
SARS-CoV were carried out in a biological safety level 3 containment labora- 


tory. 

Green fluorescence analysis. Fluorescence microscopy was used to observe the 
green fluorescence of GFP expressed from the SARS-CoV replicon. The cells 
were observed under an Olympus IX70 inverted fluorescence microscope, and 
the images were recorded using Image-Pro Plus (Media Cybernetics). 

Real-time PCR analysis. Real-time reverse transcription (RT)-PCR analyses 
were performed to quantify the copy number of SARS-CoV RNA in SCR-1 cells 
or Vero cells. Primers and RT-PCR conditions were used as previously described 
(19). Real-time PCR signals were analyzed using the LightCycler software (ver- 
sion 5.32; Roche), and the sizes and uniqueness of PCR products were verified 
by performing both melting curves and agarose gel electrophoresis. 

Statistical analysis. All data are expressed as means + standard deviations. 
Statistical significance of the in vitro data was determined by Student’s ¢ test (two 
tailed), while the significance of the differences between the median values of the 
in vivo data was determined using the two-tailed Mann-Whitney test. Statistical 
significance was assigned if P was <0.05. 


RESULTS 


Quantitative proteome analysis of SARS-CoV replicon cells. 
To gain insights into the molecular pathways perturbed by 
SARS-CoV replication, we employed SILAC combined with 
LC-MS/MS to determine the differential proteomes of the 
SARS-CoV replicon cell SCR-1 and parental BHK21 cells. 
The workflow in this study is outlined in Fig. 1. To obtain 
reliable results, we carried out both forward and reverse 
SILAC experiments. After cell lysis, SDS-PAGE fractionation, 
in-gel digestion, LC-MS/MS analysis, and the database search, 
we were able to identify 1,801 proteins, among which 1,480 
could be quantified. Among the quantified proteins, we include 
here only the quantification results for those proteins that 
could be quantified in both the forward and reverse SILAC. 
Together, this gives quantifiable results for 1,081 proteins 
(data not shown). Differentially expressed proteins (DEPs) 
were selected based on a predefined threshold of 2.0-fold 
change (the ratio of BHK21 to SCR-1 expression was greater 
than 2.0 or less than 0.5). Among the quantified proteins, 43 
were upregulated and 31 were downregulated in SCR-1 cells. 
The quantification results for the proteins with significant 
changes are summarized in Table 1. (Detailed information 
about the quantified proteins and ratios for each experiments 
was obtained [data not shown].) 

Functional categories of DEPs. In order to understand the 
biological relevance of the changes in protein expression in 
response to SARS-CoV replication, PANTHER classification 
system was used to classify the 74 DEPs according to their 
functions. The PANTHER classification system revealed that 
the DEPs can be classified into 23 groups according to their 
functional properties (Fig. 2). These proteins are implicated in 
a broad range of cellular activities (Fig. 2). Proteins involved in 
nucleic acid binding account for the largest portion (27%). 
There are also a significant number of proteins involved in 
regulatory molecule (13%), chaperone (10%), and kinase 
(9%) activities. Among these DEPs, of particular interest are 
upregulated chaperones and downregulated host translational 
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TABLE 1. List of proteins with at least 2.0-fold quantitative alteration in SARS-CoV replicon cells based on SILAC analysis 


Protein type and gene 
accession no. 


Gene product 
name 


Description 


SILAC BHK21/SCR-1 
expression ratio 


Nucleic acid binding 
IP100310317 
IPI00776411 
IPI00317794 
IPI00138335 
IP100623284 
IP100228457 
IPI00113241 
IPI00466069 
IPI00268802 
IPI00112448 
IP100474446 
IP100404707 
IPI00555113 
IP100322422 
IP100274407 
IPI00133503 
IP100756424 
IP100230679 
IP100387566 
IP100400432 


Select regulatory molecule 
IPI00131870 
IPI00127408 
IPI00315100 
IPI00754880 
IPI00114560 
IPI00762919 
IPI00113849 
IPI00138143 


Chaperone 
IP100229080 
IPI00323357 
IPI00116498 


IP100669474 
IPI00331556 
IP100230707 


IPI00331334 
IP100153740 


Kinase 
IPI00316677 
TPI00119772 
IPI00229884 
IPI00283633 
TPI00556823 
TPI00119663 


Transferase 
TPI00119772 
TPI00269091 
IPI00228457 
IPI00323353 


Transcription factor 
IPI00310317 
IPIO0119059 
IPI00135883 


Transporter 
IPI00130924 
IP100120572 
IPI00122648 


Atf3 
Hefel 
Nel 
Zfp36l1 
Sf3b1 
Myst2 
Rps19 
Eef2 
Rps18 
Rps10 
Eif2s1 
Rbm14 
Rpll8 
Mrpl47 
Tufm 
Rpl27a 
Eif5b 
Rpl36 
Ubd 
Eif4a2 


Cops3 
Racl 
Rhoa 
Arhgef12 
Rab1l 
Hip1 
Cdc42 
Cdkn2b 


Hsp90ab1 
Hspa8 
Ywhaz 


Hspala 
Hspa4 
Ywhag 


Bag3 
Ahsal 


Pik3cd 
Stk4 
Pak2 
Limk2 
Prkaal 
Mapk1 


Stk4 
Trit1 
Myst2 
Uppl 


Atf3 
Nr2fl 
Gata3 


Slc27a2 
Nup50 
Laptm4a 


Activating transcription factor 3 

Host cell factor Cl 

Nucleolin 

Zinc finger protein 36, C3H type-like 1 
Splicing factor 3b, subunit 1 

MYST histone acetyltransferase 2 

Ribosomal protein $19 

Eukaryotic translation elongation factor 2 
Ribosomal protein $18 

Ribosomal protein $10 

Eukaryotic translation initiation factor 2, subunit 1 alpha 
RNA binding motif protein 14 

Ribosomal protein L18 

Mitochondrial ribosomal protein L47 

Tu translation elongation factor, mitochondrial 
Ribosomal protein L27a 

Eukaryotic translation initiation factor 5B 
Ribosomal protein L36 

Ubiquitin D 

Eukaryotic translation initiation factor 4A2 


COP9 signalosome complex subunit 3 
RAS-related C3 botulinum substrate 1 

Ras homolog gene family, member A 

Rho guanine nucleotide exchange factor (GEF) 12 
RAB1, member RAS oncogene family 

Huntingtin interacting protein 1 

Cell division cycle 42 

Cyclin-dependent kinase inhibitor 2B 


Heat shock protein 90 kDa alpha (cytosolic) 

Heat shock protein 8 

Tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation 
protein, zeta polypeptide 

Heat shock protein 1A 

Heat shock protein 4 

3-Monooxygenase/tryptophan 5-monooxygenase activation protein, 
gamma polypeptide 

Bcl2-associated athanogene 3 

AHAI, activator of heat shock protein ATPase homolog 1 (yeast) 


Phosphatidylinositol 3-kinase catalytic delta polypeptide 
Serine/threonine kinase 4 

P21 (CDKN1A)-activated kinase 2 

LIM motif-containing protein kinase 2 

Protein kinase, AMP-activated, alpha 1 catalytic subunit 
Mitogen-activated protein kinase 1 


Serine/threonine kinase 4 

TRNA isopentenyltransferase 1 
MYST histone acetyltransferase 2 
Uridine phosphorylase 1 


Activating transcription factor 3 
Nuclear receptor subfamily 2, group F, member 1 
GATA binding protein 3 


Solute carrier family 27 (fatty acid transporter), member 2 
Nucleoporin 50 
Lysosomal-associated protein transmembrane 4A 


0.18 + 0.08 
0.19 + 0.04 
0.25 + 0.06 
0.29 + 0.05 
0.30 + 0.03 
0.41 + 0.04 
2.27 + 0.27 
2.29 + 0.34 
2.35 + 0.41 
2.46 + 0.71 
2.83 + 0.53 
2.96 + 0.25 
3.07 + 0.19 
3.16 + 0.39 
3.18 + 0.41 
3.51 + 0.77 
4.08 + 0.60 
4.08 + 0.71 
5.40 + 0.85 
733-2 1.77 


0.21 + 0.09 
0.22 + 0.09 
0.31 + 0.05 
0.40 + 0.09 
0.46 + 0.08 
3.50 + 0.52 
4.16 + 0.15 
4.58 + 0.70 


0.11 + 0.03 
0.18 + 0.05 
0.23 + 0.04 


0.24 + 0.07 
0.31 + 0.06 
0.31 + 0.07 


0.37 + 0.04 
0.41 + 0.07 


0.18 + 0.01 
0.24 + 0.06 
0.25 + 0.03 
0.25 + 0.06 
0.29 + 0.05 
0.29 + 0.07 


0.24 + 0.06 
0.41 + 0.03 
0.41 + 0.04 
2:37 £0.37 


0.18 + 0.08 
0.31 + 0.10 
0.32 + 0.04 


0.33 + 0.05 
0.33 + 0.15 
4.53 + 1.03 


Continued on following page 
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Protein type and gene 
accession no. 


Gene product 
name 


Description 


SILAC BHK21/SCR-1 


expression ratio 


Cytoskeletal 

IPI00515564 Kifcl Kinesin family member Cl 2.90 + 0.34 

IPI00762919 Hip1 Huntingtin interacting protein 1 3.50 + 0.53 

IPI00229647 Tin2 Talin 2 5.33 + 0.55 
Oxidoreductase 

IPI00318108 Acox3 Acyl-coenzyme A oxidase 3, pristanoyl 2.42 + 0.44 

IP100229510 Ldhb Lactate dehydrogenase B 2.43 + 0.37 

IPI00222809 Hopd Hexose-6-phosphate dehydrogenase (glucose 1-dehydrogenase) 4.77 + 0.49 
Signaling molecule 

IP100320634 Elk3 ELK3, member of ETS oncogene family 0.20 + 0.05 
Synthase and synthetase 

IP100331707 Hmgcs1 3-Hydroxy-3-methylglutaryl-coenzyme A synthase 1 0.21 + 0.03 
Select calcium binding 

protein 

IP100317309 Anxa5 Annexin A5 0.43 + 0.09 
Receptor 

IPI00119059 Nr2f1 Nuclear receptor subfamily 2, group F, member 1 0.30 + 0.10 
Membrane trafficking 

IP100266752 Cpne3 Copine-3 0.22 + 0.06 
Lyase 

IPI00318496 Gad1 Glutamic acid decarboxylase 1 0.30 + 0.03 

IP100462072 Enol Enolase 1, alpha non-neuron 2.06 + 0.13 
Protease 

IPI00128154 Ctsl Cathepsin L 2.78 + 0.46 

IP100626909 Capn1 Calpain 1 4.32 + 0.88 
Phosphatase 

IP100116554 Ptpn11 Protein tyrosine phosphatase, non-receptor type 11 0.28 + 0.04 

IPI00130507 Dusp7 Dual-specificity phosphatase 7 2.58 + 0.35 

IPI00330483 Ppp1r10 Protein phosphatase 1, regulatory subunit 10 3.06 + 0.39 

IP100177072 Pfkfb3 6-Phosphofructo-2-kinase/fructose-2,6-biphosphatase 3 3.09 + 0.33 
Miscellaneous function 

IPI00222515 Psmd11 Proteasome (prosome, macropain) 26S subunit, non-ATPase, 11 0.26 + 0.05 

IP100669522 Clasp1 CLIP associating protein 1 0.34 + 0.06 
Molecular function 

unclassified 

IPI00169634 Snapcl Small nuclear RNA activating complex, polypeptide 1 0.21 + 0.08 

IP100315187 2400001 E08Rik RIKEN cDNA 2400001E08 gene 0.22 + 0.03 

IPI00126917 Phb Prohibitin 0.37 + 0.06 

IP100654076 F630110N24Rik RIKEN cDNA F630110N24 gene 0.40 + 0.07 

IPI00187289 2310028N02Rik RIKEN cDNA 2310028N02 gene 3.00 + 0.23 


machinery. We expect that these host factors play important 
roles in SARS-CoV replication. 

Validation of differential protein expression. To further con- 
firm the SILAC ratios we observed with MS, we used Western 
blot analysis to examine the expression of a selected panel of 
proteins, including HSPA1A, HSP90AB1, hnRNPAI, lactate 
dehydrogenase B (LDHB), RPS19, YWHAZ, ENO1, cell di- 
vision control protein 42 (CDC42), and BAG3. As shown in 
Fig. 3, for all of the selected proteins, the Western blotting and 
densitometric analysis results showed the same pattern of ex- 
pression as that obtained from SILAC experiments. 

BAG3 is essential for efficient SARS-CoV replication. To 
evaluate the functional requirement of BAG3 for the replica- 


tion of SARS-CoV, we used RNA interference (RNAi) to 
reduce cellular BAG3 levels and examine the effect of BAG3 
knockdown. The effect of the BAG3 siRNA in silencing the 
BAG3 gene in SCR-1 cells was examined by directly measuring 
changes in BAG3 protein levels. As shown in Fig. 4B, trans- 
fection of replicon cells with the BAG3 siRNA resulted in a 
significant decrease in the level of BAG3 protein, whereas 
mock transfection with a random siRNA had no such effect. 
This effect was specific because the BAG3 siRNA did not 
change the levels of tubulin protein (Fig. 4A). 

To analyze whether BAG3 knockdown has an effect on 
SARS-CoV replication in SCR-1, the green fluorescence anal- 
ysis, quantitative real-time PCR, and Western blotting were 
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FIG. 2. Pie chart representations of the distribution of differentially 
expressed proteins according to their molecular functions (A) and 
biological processes (B). Categorizations were based on information 
provided by the online resource the PANTHER classification system. 


performed. As shown in Fig. 4, knockdown of BAG3 resulted 
in a marked decrease in GFP fluorescence signal, replicon 
RNA levels, and SARS-CoV nsp5 protein level in SCR-1 cells. 
All of these data indicate that BAG3 knockdown can inhibit 
SARS-CoV replication and protein synthesis. 

To evaluate the functional requirement of BAG3 for the 
replication of SARS-CoV in Vero E6 cells, we constructed cell 
lines stably expressing siRNAs targeting the BAG3 mRNA. 
The knockdown of BAG3 in Vero-KD cells was confirmed by 
Western blotting analysis (Fig. 5A). The apoptosis and cell 
growth assay revealed that reduction of BAG3 had no signif- 
icant influence on Vero cell growth or apoptosis (data not 
shown). The resulting cell lines were tested for the ability of 
SARS-CoV to replicate using real-time PCR and plaque re- 
duction assay. Reduction of BAG3 levels resulted in a signif- 
icant decrease in SARS-CoV RNA and virus titer compared to 
the parental Vero E6 and Vero-NC cells (Fig. 5B and C). To 
evaluate BAG3’s specificity for SARS-CoV replication, we 
asked if growth of another coronavirus, IBV, was altered in the 
Vero-KD cells. In contrast to SARS-CoV, plaquing efficiency 
of IBV was unaffected by decreased BAG3 levels (Fig. 5C). 
Thus, our data indicate that BAG3 is a host protein that is 
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A Protein SILAC Ratio Western Results B Normalized Fold Change 


Name (BHK21/SCR-1)” BHK21 SCR-1 4 2 #0 2 4 
HSP90AB1 0.11 —-= HSP90AB1 
HNRNPAI 0.24 -= HNRNPA1 
HSPA1A 0.24 —— HSPAIA 
LDHB 2.43 -— LDHB 
RPS19 2.28 _—— RPS19 
cDC42 4.16 o--—| cDC42 
ENO1 2.06 —— ENO1 
YWHAZ 0.24 | YWHAZ 
BAG3 0.37 —— BAG3 
Tubulin 1.04 oe 


FIG. 3. Western blot and densitometric analysis of nine DEPs and 
internal control protein tubulin. (A) Western blot images for nine 
DEPs and internal control protein tubulin. (B) Relative expression of 
nine DEPs (normalized to tubulin band) was determined using image 
densitometry and expressed in a bar chart format. The error bars 
represent the standard deviation of the mean for three representative 
analyses. 


specifically required for efficient replication of SARS-CoV and 
not IBV. 


DISCUSSION 


The emergence of SARS-CoV underscores the importance 
of advancing understanding of virus-host interactions. Increas- 
ing evidence emphasizes comparative proteomics to screen the 
differentially expressed proteins associated with host cellular 
pathophysiological processes of virus infection (42, 54, 63). 
From the literature, very few studies have been performed to 
analyze the interaction between CoV and host cells using pro- 
teomics analysis. This study is the first to employ the SILAC 
technique to globally search for the dysregulated host proteins 
in SARS-CoV replicon cells. In the present work, a total of 74 
differentially expressed host proteins were identified and dif- 
ferential expression levels of nine DEPs were confirmed by 
Western blotting and densitometric analysis. The PANTHER 
classification system revealed that the proteins can be classified 
into 23 groups according to their biological process or molec- 
ular functions (Fig. 2). Based on the identified proteins in the 
present work, we obtained an overview of the altered protein 
expression of host cells responding to SARS-CoV replication. 

By comparison with the previous genomic study that ana- 
lyzed the SARS-CoV-infected cells using mRNA microarray 
approaches (37), among 74 proteins identified in the current 
study, 14 of them were also found to be altered at the mRNA 
level, namely, Rps10, Rps19, Laptm4a, Racl, Ctsl, Rpl27a, 
Myst2, Rpl36, Rps18, Tritl, Ppp1r10, Rpll8, Hspala, and 
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FIG. 4. Effects of BAG3 knockdown on SARS-CoV replicon cells. 
(A) The knockdown of BAG3 in SCR-1 cells was confirmed by West- 
ern blotting analysis. Western blotting results indicated that BAG3 is 
upregulated in SCR-1 cells, and knockdown of BAG3 led to reduction 
of SARS-CoV proteins in SCR-1 cells. (B) Inhibition of BAG3 expres- 
sion led to suppressed GFP fluorescence in SCR-1 cells compared to 
untreated SCR-1 cells or cells transfected with control siRNA. 
(C) Quantification of SARS-CoV mRNA revealed markedly reduced 
mRNA in BAG3 knockdown SCR-1 cells compared to untreated cells 
or cells transfected with control siRNA. 
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Clasp1. Most of the remaining identified proteins are newly 
discovered. The limited overlap between two studies is ex- 
pected to reflect, at least in part, the different effects between 
viral infection and replicon RNA on the host cell. In addition, 
the poor correlation between proteomic and genomic results 
was not unexpected, as previous studies showed that there 
exists an approximately 60 to 80% discordance between 
mRNA and protein abundances (47, 64). This lack of correla- 
tions could be the result of mRNA degradation, alternative 
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FIG. 5. BAG3 is essential for efficient SARS-CoV replication. Col- 
umns are the means of three independent experiments; error bars 
represent + standard deviation. The asterisks indicate significant dif- 
ference (P < 0.05). (A) Confirmation of BAG3 knockdown in 
Vero-KD cells by Western blotting analysis. (B) Quantification of 
SARS-CoV RNA revealed markedly reduced SARS-CoV RNA in 
Vero-KD cells compared to parental Vero E6 or Vero-NC cells. Quan- 
titative PCR was performed as described in Materials and Methods. 
(C) The reduction in SARS-CoV or IBV titers was determined by 
plaque assays. Vero E6, Vero-NC, and Vero-KD cells were infected 
with SARS-CoV or IBV, and the plaque numbers in the different cell 
lines were then normalized to the number of plaques that formed in 
the Vero E6 cells. 


splicing, and posttranscriptional regulation of gene expression. 
Therefore, the current proteomic study represents a comple- 
mentary strategy for studies of mechanisms that underlie the 
replication of SARS-CoV. 

Among all DEPs, significantly upregulated BAG3, a mem- 
ber of the BAG cochaperone protein family, was selected for 
further functional studies. Our results demonstrate that BAG3 
knockdown results in greatly diminished virus replication, sug- 
gesting the possibility that upregulation of BAG3 may be part 
of the machinery that SARS-CoV relies on for replication. 
However, the mechanisms by which BAG3 regulates SARS- 
CoV replication remain unknown. BAG proteins participate in 
a wide variety of cellular processes, including cell survival, 
cellular stress response, proliferation, migration, and apoptosis 
(14, 25). BAG3 is implicated in the pathogenesis of neoplasia 
via its ability to regulate stress-induced apoptosis in a prosur- 
vival fashion. This regulation occurs at a number of levels, 
including cytochrome c release, apoptosome assembly, and 
others (5, 55). Other cellular signaling molecules that have 
been reported to be regulated by BAG3 include Raf-1, cyclin- 
dependent kinase 4 (CDK-4), and epidermal growth factor 
receptor (EGFR) (13), as well as focal adhesion kinase (27). 
Notably, there is increasing evidence that BAG3 is required for 
efficient growth of different viruses, including varicella-zoster 
virus (VZV) (35), polyomavirus JC (3), Epstein-Barr virus 
(EBV) (66), herpes simplex virus (HSV) (36), and HIV (56) 
and depletion of BAG3 by small interfering RNA results in 
inhibition of virus replication (35). 
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Recent evidence implicates an additional function of BAG3 
in the regulation of the autophagy pathway. These findings 
indicate that autophagosome formation and turnover may de- 
pend on BAG3 and BAG3 can stimulate autophagy processes 
(7, 8, 18). Autophagy is a cellular stress response that functions 
to recycle proteins and organelles and an intracellular cata- 
bolic transport route conserved among all eukaryotic cells (22, 
65). One of its functions is to act as an immune mechanism 
against intracellular pathogens (38, 50). Bacteria and viruses 
targeted for destruction are sequestered into large double- 
membrane vesicles (DMV) called autophagosomes and subse- 
quently delivered to the lysosomes, where they are degraded by 
hydrolases. However, it has become clear that the autophagy 
pathway is often exploited by viruses to facilitate their entry or 
replication (29, 32). Importantly, the laboratory of Mark Deni- 
son has shown that the DMVs induced by infection with 
SARS-CoV and the coronavirus mouse hepatitis virus (MHV) 
could be autophagosomes, while MHV infection appears to 
induce autophagy and inhibition of autophagy leads to sup- 
pression of MHV growth (52, 53, 68). Their results demon- 
strate that the autophagy pathway is required for formation of 
a viral replication complex and for efficient viral growth. The 
significant conservation of replicase gene organization and 
probable functions of SARS-CoV and MHV replicase proteins 
further suggest that SARS-CoV may share with MHV the 
strategy of induction and use of the autophagy pathway for 
efficient replication (53). 

Based on the critical role of BAG3 in the stimulation of the 
autophagy pathway, it is tempting to suggest that BAG3 might 
be recruited to and activate the autophagy machinery to facil- 
itate the efficient replication of SARS-CoV. Therefore, the 
increase in BAG3 expression is likely part of the cell’s response 
to SARS-CoV infection and appears to represent a novel 
mechanism for maintaining SARS-CoV replication in host 
cells. This speculative idea, however, is not yet supported by 
experimental data, and further experiments will be required to 
determine the functional implication of BAG3 and the auto- 
phagy pathway in the life cycle of SARS-CoV. 

In the present work, the most distinguishable proteins to be 
downregulated in SARS-CoV replicon cells were those in- 
volved in the host translational machinery, including 40S ribo- 
somal proteins (RPS10, RPS18, and RPS19), 60S ribosomal 
proteins (RPL18, RPL27A, and RPL36), and eukaryotic trans- 
lation initiation factors (EIF2S1, EIF5B, and EIF4A2). Such 
host translational downregulation is frequently observed, as 
evident in infections with herpes simplex virus (15), poliovirus 
(34), and SARS-CoV (37). Recently, it was reported that 
SARS-CoV nsp1 protein uses a two-pronged strategy to inhibit 
host translation and gene expression by binding to the 40S 
ribosomal subunit and inactivating the translational activity of 
the 40S subunits (26). Thus, downregulation of the host trans- 
lational machinery may play an important role in SARS-CoV 
replication in host cells. 

Two regulatory molecules, Cdc42 and RhoA, were found to 
be downregulated in replicon cells. Cdc42 and RhoA belong to 
the Rho family of small GTPases. These proteins play a role in 
cell migration, membrane traffic, and actin cytoskeleton reor- 
ganization (6). It has also been shown that Cdc42 regulates the 
phosphatidylinositol 3-kinase (PI3K)-mTOR pathway (6, 16) 
and PI3K/Akt signaling pathways are important for establish- 
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ing persistent SARS-CoV infection in Vero E6 cells (44, 45). 
Therefore, modulation of SARS-CoV replication by Cdc42 
and RhoA may be mediated by the PI3K-mTOR pathway and 
the PI3K-mTOR pathway may be one of the key factors for 
understanding persistence of SARS-CoV replicon RNA in 
host cells. 

Other proteins of interest uncovered in this study include 
14-3-3 proteins. The 14-3-3 proteins have emerged as critical 
regulators of diverse cellular responses in eukaryotic organ- 
isms (2, 17, 46). In mammalian cells, seven different isoforms 
have been identified (¢, B, y, €, 0, n, and 8), with each isoform 
having distinct tissue localization and function. Through inter- 
actions with more than 400 target proteins identified so far, 
these 14-3-3 proteins are known to be involved in widespread 
biological processes such as signal transduction, cell cycle con- 
trol, apoptosis, cellular metabolism, proliferation, cytoskeletal 
regulation, transcription, redox regulation, stress response, etc. 
(11, 61). Interesting evidence has been reported indicating that 
14-3-3 might be important for host translational shutoff, the 
shutoff of minus strand synthesis, or other processes that are 
time dependent in virus-infected cells (12). Recently, it was 
also reported that 14-3-3 proteins are involved in the activation 
of signaling cascades during viral infection (48). Therefore, the 
increase in 14-3-3 protein expression is likely part of the cell’s 
response to SARS-CoV replication and appears to represent a 
homeostatic mechanism for cell defense. 

In conclusion, the quantitative proteomics analysis described 
here identified many host factors that potentially affect SARS- 
CoV replication and implicated previously unconsidered path- 
ways in the virus replication. By using siRNA, we also obtained 
evidence that depletion of BAG3 results in inhibition of virus 
replication. Our study provides important information for un- 
derstanding SARS-CoV replication. Further studies will deter- 
mine more directly how implicated host factors affect the virus 
and how such effects illuminate cellular functions and path- 
ways. 
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